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ABSTRACT 

Evidence is presented indicating that the bursting component of the X-ray radiation 
detected in the nuclear region of the active radio galaxy 3C 111 comes from the blobs 
ejected in the pc-scale jet and not from the accretion disc. After each new outburst 
the radio flux density associated with it increases to a peak in ~ 1 year and then 
subsides over a period of 1-2 years with the flux falling off exponentially as the blob 
moves outward and dissipates. Similar peaks (bursts) are seen in the X-ray light curve 
and a cross-correlation between the two shows a very high correlation with the X-ray 
peaks leading the radio peaks by ~ 100 days. A second cross-correlation, this time 
between the radio event start times and the X-ray light curve, also shows a significant 
correlation. When this is taken together with the long (~ 1 yr) delay between the 
start of each ejection event and its associated X-ray peak it indicates that this bursting 
component of the X-ray flux must be associated with the ejected blobs in the pc-scale 
jet and not with the accretion disc. Because X-ray telescopes do not have the resolution 
required to resolve the accretion disc area from the pc-scale jet, this paper is the first 
to present observational evidence that can pinpoint the point of origin of at least those 
long-timescale X-ray bursts with durations of 1-3 yrs. 

Subject headings: galaxies: active — galaxies: jets — X-rays: bursts 



1. Introduction 

Variable radio loud active galactic nuclei 
(AGN) with their associated jets have been 
studied for more than 40 years. At radio fre- 
quencies tremendous advances in their under- 
standing came with the better than 1 mil- 
liarcsec resolution of the VLBA which al- 
lowed disturbances propelled outward from 
near the central black hole to be tracked as 
they moved rapidly outward in a tightly con- 
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fined jet. Similarly, advances came at X-rays 
when the Chandra telescope came on-line. 
Unfortunately, the resolution at X-ray wave- 
lengths is much lower (0.5 arcsec) than with 
the VLBA. After 40 years of study the ex- 
act mechanism by which the jets are formed 
and the ejected material is accelerated is still 
unclear. Even the exact point of origin of 
the non-varying radio component of the AGN 
core rad iation is still bei n g de bated, with 
for e.g. Bell and Comeau ( 20ld ) presenting 
arguments favoring a core component cen- 
tered on the accretio n disc and others e.g. 
Marscher etHI (|2010h arguing for . 



a core com- 



ponent located outside the accretion disc area 
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in the base of the jet. 

Variations in the total flux of AGN sources 



at rad io frequencies were studied by lAndrew et al. 



( 19781 ) when large, single-dish telescopes were 
first coming on line. Before it was possible to 
map these sources with the milliarcsec reso- 
lution of th e VLBA it was demonstrated by 



Lega (|1984 see Fig 1 below) that the flux 



bursts from different AGN sources, though 
differing slightly in time scale, had similarly 
shaped profiles. In each event the flux was 
observed to increase to a peak within ~ 1 
year and then decrease exponentially, falling 
below the detection limit within 1-3 years. 
Now, using the VLBA, individual bursts can 
be tracked at radio frequencies as separate 
plasma blobs that move away from the cen- 
tral object, reach a pe ak and then slowly dis- 



sipate over 1-3 years (jChatterjee et al.l 12011 
see their Fig 6). 

Before a complete understanding of what 
is happening in the nuclear region of these 
sources is likely to be achieved it is essen- 
tial that the exact point of origin of the dif- 
ferent radiation components be determined. 
One of the biggest uncertainties remaining is 
where the bursting component of_the X-ra 



radiation comes from in AGNs (jGallol [2006; 
Beckmann et al. 201 ll ). Is it from material 



ass ociated w i th the accretion disc, as argued 
bv lMarscherl (l2006h . or is it from the pc-scale 
jet? These components cannot be resolved 
by Chandra. Although some X-ray radia- 
tion has been detected in the ou ter, kpc-scale 



jet an d terminal knot of 3C 111 (jHogan et al 



l201ll ). most comes from the unresolved core 
region containing the accretion disc and pc- 
scale jet. Much work has been done examining 
Seyferts and X-ray binaries (IZdziarski et a] 



2003; iGallo et alJbooi iKing et al 



2011 



but 

in most cases these observations involve vari- 
ations with shorter timescales (less than 0.3 
yr) than those being examined here. Be- 
fore proceeding, it is important to define ex- 



actly what timescale is referred to here when 
the term "burst" is used. In our case the 
burst timescale is that defined in early vari- 
able source work and covers a period of 1- 
3 yrs, as defined by a typical radio outburst 
at cm wavelengths. This has been shown by 
VLBA observations to be the time taken for 
the radiation from a new radio outburst to 
reach a peak and dissipate as it moves out in 
the pc-scale jet. There are many bursts at 
all wavelengths that have shorter time scales 
but this likely means that they are occurring 
in relatively small structures near the super- 
massive black hole. 



From Fig 14 of lChatteriee et al.l (|201ll ) the 
X-ray radiation from 3C 111 can be seen 
to be composed of two main components of 
roughly equal peak strengths. These are: 1), 
a steady component and 2), a superimposed 
bursting component. The former is assumed 
to be associated with the accretion disc area. 
It is this component that may be quenched 
briefly just prior to the start of each new radio 
event, producing the di p that was reported by 
Chatteriee et all (|201lh . Its appearance close 



to the event start time requires that it orig- 
inate close to the nucleus. The latter vari- 
able component, is composed of several X-ray 
bursts that are similar to, or slightly brighter 
than, the steady component, and it may be 
possible to determine the point of origin of 
this component by carrying out cross corre- 
lations a) between the radio and X-ray light 
curves and b) between the radio event start 
times and the peaks in the X-ray light curve. 
Although there may be a steady X-ray compo- 
nent associated with the accretion disc, here 
we are concerned mainly with the bursting 
one, and we show in Sections 3 and 4 that 
most of the bursting component of the X-ray 
continuum in the FR-II radio galaxy 3C 111 
originates in the plasma blobs in the pc-scale 
jet as opposed to the accretion disc area. 
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2. The Data 
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Fig. 1. — Typical radio burst profile from 
Legg (1984). 



Recently I Chatter jee et al.l (120111 ) have pre- 
sented a significant study of the bursting com- 
ponents of the radiation from 3C 111 at sev- 
eral frequencies ranging from 15 GHz to X- 
rays. This included VLBA maps of the jet 
structure at 37 GHz for ~ 50 different epochs 
between 2004 and 2010. At radio frequencies, 
with the milliarcsec resolution provided by the 
VLBA the progression of each newly ejected 
blob is easily followed as it moves outward in 
the inner pc-scale jet. Where it was only as- 
sumed in the early days that the variable com- 
ponent of the radio flux was associated with 
the ejected blobs, it can now easily be seen 
from a series of consecutive VLBA maps that 
this is the case. Also, 3C 111 is a much better 
source to use for this purpose than is 3C 120, 
BL Lac, or PKS 1510-089 whose outbursts are 
so frequent that their light curves are much 
more confused. 



Chatteriee et al.l (|201ll ) also reported see 



ing a temporary dip in the X-ray radiation 
just prior to the ejection of each new radio 
event. Although not a perfect correlation 
these dips do appear convincing and may be 
due to a decrease in the otherwise non- varying 
X-ray radiation originating near the accretion 
disc. Although the timing of these dips pre- 
cedes the estimated event start times by ~ 50 
days, there is no mention of whether any al- 
lowance was made for an acceleration period 
between the actual event time and the first ap- 
pearance of a new blob. In fact, the dips may 
occur much closer to the actual event start 
time. It is assumed here that they may be due 
to some kind of initial quenching that occurs 
in the non-varying X-ray component associ- 
ated with the accretion disc when new mate- 
rial is first injected. The presence of an X-ray 
dip appearing shortly prior t o the ejection of 



a radi o blob was also used bv lMarscher et al 



([20021 ) ;o argue for a connection between the 
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black hole and the jet in 3C 120. Unfortu- 
nately, as noted above, because of the poor 
resolution at X-rays, it is not possible to pin- 
point the point of origin of the X-ray ra- 
diation. We can assume, however, that X- 
ray variations that correlate closely in time 
(within a few tens of days) with the start of 
ejection events are likely to be tied to the ac- 
cretion disc area while peaks in the X-ray light 
curve that correlate with radio ejection events, 
but appear a long time after the start of the 
radio event itself, must originate much further 
out in the jet. 

In Fig 2 the solid curve rep resents the X-ray 



data o btained for 3C 111 bv lChatteriee et al 



and was reproduced here from a dig- 
itized version of their smoothed curve. The 
times at which each new radio event began 
are indicated by a vertical bar and the let- 
ter E. Also in Fig 2 the dips reported by 
Chatteriee et al.1 (l201lh have been highlighted 
and the dashed curves give a rough indication 
of how the light curve would be expected to 
look if the dips were not present. 

When the light curves of 3C 111 ob- 
tained at differen t wave lengths are compared 
(jChatteriee et al.l 1201 ll . see their Fig 4), it 



is apparent that for each burst the shorter 
wavelengths lead those at radio, with the X- 
ray radiation leading the radio by ~ 100 days. 
Using this information together with the pro- 
file shape from Fig 1, if the X-ray radiation 
comes from the jet, we can predict approxi- 
mately when each X-ray peak will be expected 
to occur, at least for each peak that occurs 
after 2005.0, and these times are indicated in 
Fig 2 by P. It is immediately apparent that 
there is good agreement between the predicted 
and observed X-ray peaks for all cases except 
during 2008. Here, when E6 is predicted to 
be peaking, three new events (E7, E8, and 
E9) are about to begin and the previously 
reported dips may also be occurring. In ad- 
dition to this, for a long period (2008.3 to 




Fig. 2. — Plot of X-ray flux between 2004 and 
2010 from Chatterjee et al.(2011). Dips dis- 
cussed by Chatterjee are shown shaded and 
are indicated by D. Times when new events 
occurred are indicated by E while the times 
that the flux from these events is estimated 
to peak, from Fig 1 of Legg (1984), are indi- 
cated by P. Estimated burst profiles are indi- 
cated by the dashed curves. Events 2 and 3, 4 
and 5, and 7, 8, and 9 are considered as single 
events for the purpose of defining profiles. 
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2008.6) no data are available. 

There appears to be little evidence that in 
each burst the increase in flux and its drop- 
off again, as outlined by the dashed curves, is 
related to Doppler boosting, since each blob 
appears to move out with a consta nt velocity 
in Fig 7 of IChatteriee et all (|201lh . 



3. Cross-Correlation Between 14.5 GHz 
Radio and X-ray Light Curves 

In Fig 3 the dips highlighted in Fig 2 have 
been removed. Since they appear close in time 
to the new event, they must originate near the 
accretion disc and will not then be part of the 
jet radiation. Since their presence could con- 
fuse the cross-correlation analysis between the 
event times and the light curves from the jet 
that we want to examine, they have been re- 
moved prior to the initial calculation. Note 
that a second analysis using data that in- 
cluded the dips was also carried out to insure 
that the removal of the dips was not signifi- 
cantly affecting the correlation results. 

Since the bursting component at radio fre- 
quencies is known, from the VLBA maps, to 
come from the blobs ejected in the pc-scale 
jet, a strong correlation between the radio 
(14.5 GHz) and X-ray light curves would be 
strong evidence that the bursting X-ray radi- 
ation also originates in the pc-scale jet. We 
therefore carried out a cross-correlation anal- 
ys is between these t wo cu rves shown in Fig 4 
of lChatteriee et all (1201 lh . To do this the two 
curves were first digitized into bins of length 
0.05yr. The correlation coefficient was then 
calculated between the two curves as one was 
shifted relative to the other. 

In Fig 4 the squared Pearson correlation 
coefficient obtained in this analysis is plotted 
as a function of the delay in days between the 
radio and X-ray peaks. The solid curve was 
obtained using only that data in Fig 3 before 
yr 2008.3. This portion of the data was looked 
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Fig. 3. — Same as Fig 2 but with dips re- 
moved. 
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Fig. 4. — Squared correlation coefficient be- 
tween X-ray and 14.5 GHz radio light curves 
plotted as a function of delay in number of 
days that radio follows the X-ray. (solid 
curve) data prior to yr 2008.3. (dashed curve) 
all data. 



at first because it is difficult to know how the 
previously mentioned 4-month data gap in yr 
2008, and the presence of three consecutive 
dips (associated with radio events E7, E8, and 
E9) would affect the results. For the data be- 
fore 2008.3 a very significant correlation is ob- 
tained (P = 2 x 10 -4 ). The result also indi- 
cates that the radio peaks are delayed by ap- 
proximately 100 days. The squared Pearson 
correlation coefficient obtained using all the 
data in Fig 3 is shown by the dashed curve 
in Fig 4. Only a slightly poorer probability 
(P = 5 x 10 -4 ) is obtained. The correlation 
analysis was then repeated using the data in 
Fig 2 which contains the dips. The result 
was only marginally worse and still resulted 
in a highly significant probability. This in- 
dicates that the strong correlation found has 
not been produced by the removal of the dips 
but is, instead, only slightly strengthened, as 
predicted. This further indicates that the ac- 
curacy of the procedure used to remove the 
dips cannot significantly affect the correlation 
results. This gives us strong evidence to argue 
that the radio and X-ray peaks are truly cor- 
related, which in turn is grounds to argue that 
the radio event start times may also be corre- 
lated with the X-ray peaks if the appropriate 
time delay is used. 

4. Cross-Correlation Between Radio 
Event Start Times and Peaks in the 
X-ray Light Curve 

We then carried out a cross-correlation 
analysis between the start times of the radio 
events and the X-ray light curve. The results 
are shown by the solid curve in Fig 5. Here 
the correlation coefficient has been obtained 
between the radio event start times and the 
structure in the X-ray light curve in Fig 3, as 
one curve is shifted relative to the other. As 
predicted from the curve in Fig 1, a correla- 
tion peak is seen between the radio event start 
times and the peaks in the X-ray light curve 



6 



that occurs when the radio curve is shifted 
by 328 days relative to the X-ray light curve. 
This corresponds to a shift of 17.9 channels, 
where the channel width is set by the bin size 
of 0.05 yr (18.3 days) used in the digitiza- 
tion process. The probability represented by 
this peak is P = 0.013 and is quite signifi- 
cant. But this should be no surprise since it 
is already obvious from Figs 2 and 3 that the 
predicted P-values align well with the peaks 
in the X-ray light curve. Although there is 
a weaker correlation peak obtained when the 
radio and X-ray curves are shifted relative to 
each other by channel shifts between 2 to 6 
in Fig 5, it is negative in r, where r is the 
correlation coefficient, and only becomes pos- 
itive when r is squared. This peak represents 
an anti-correlation associated with the valleys 
between the peaks. A cross-correlation was 
also carried out between the radio event start 
times and the data in Fig 2, which still con- 
tains the X-ray dips. Again it showed only a 
marginal worsening in the correlation coeffi- 
cient when the dips were included. Unlike the 
negative dip components, which must origi- 
nate in material near the accretion disc, there 
is now little doubt that the positive, bursting 
X-ray components arise in the plasma blobs 
in the pc-scale jet, well outside the accretion 
disc. 

5. Discussion 

Although there is no evidence for a corre- 
lation peak in Fig 5 near -2 channels, which 
would be produced by the dips reported by 
Chatteriee et al.l ( 201ll ). this is expected since 
the analysis has been carried out on the data 
in Fig 3 where the dips have been removed. In 
Fig 6 a cross-correlation between event times 
and the data in Fig 2, which still contains the 
dips, is given by the solid curve. There is now 
evidence for a peak, although with a probabil- 
ity of P = 0.061 it is only marginally signifi- 
cant. It is found to occur ~ 36 days prior to 
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Fig. 5. — Squared correlation coefficient for 
data in Fig 3. Event times are shifted and cor- 
related with the X-ray light curve. The peak 
between 2 and 6 channels is anti-correlated 
(negative) but inverted when the correlation 
coefficient is squared. It is associated with 
the valleys between the bursts. See text for 
an explanation of the dashed curve. 



7 



0.04 



I channel = 18 25 days 




-8 -6 -4-2 2 

Shift (channels) 

Fig. 6. — (solid curve) Squared correlation co- 
efficient for data in Fig 2. Event times are 
shifted and correlated with the X-ray light 
curve. Here the peak is negative prior to 
squaring because the dip is negative. Verti- 
cal bar indicates location of the correlation 
peak found by Chatterjee with dips 6 and 9 
omitted from the calculation, (dashed curve) 
Same for data in Fig 2 with low weight given 
to event 6. 
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Fig. 7.— Fe Ka line (solid) and F(2-10) con- 
tinuum (dashed) plotted versus date. Each 
peak can be seen to rise first in the contin- 
uum. Data are from Chatterjee (2011). 
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Fig. 8. — Pearson correlation coefficient plot- 
ted versus shift of Fe Ka line relative to F(2- 
10) continuum. The curve center indicates 
that the Fe Ka line variations lag ~50 days 
behind those of the continuum. 



the ev ent start times where IChatteriee et al 



l|201ll ) bund a value of 54 days. This may b 



due to the fact that these authors did not in- 
clude events 6 and 9 in their calculations. In 
fact, there seems to be some question about 
the reality of b u rst E 6 which, in Fig 7 of 
Chatteriee et al.1 (|201lh . appears to be con- 
fused with E5. 

To examine the effect that events 6 and 9 
might have on the correlation results, the cor- 
relation was re-run giving each of these events 
low weight in turn. It made no difference to 
the result in Fig 6 whether event 9 was in- 
cluded in the calculation or left out. How- 
ever, when event 6 was given low weight the 



significance was greatly enhanced resulting in 
a correlation coefficient of r 2 = 0.035 and a 
probability of P = 0.019, as is shown by the 
dashed curve in Fig. 6. However, it did not 
change the time found between the dip and 
its associated event. 

Giving burst E6 low weight in the cross- 
correlation calculation in Fig 5 resulted in 
only a marginal increase in significance as is 
shown by the dashed curve in Fig 5. 

X-rays have been found previously to be 
associated with radio blobs in the kpc-scale 
jet of 3C 111 and other sources where the 
blobs can b e resolved by Chandra as well as 



the V LBA (jHogan et al.ll201ll ; iKataoka et al 



20031 ). That work has established that, al- 



though the radiation from the core region al- 
ways dominates, X-ray emission from large 
scale jets of radio galaxies and quasars is more 
common than expected with X-ray jets de- 
tected in 78% of the sources studied. 

The results found here contradict what is 
usually assumed to be the case, namely that 
the bursting X-ray component is associated 



with the accretion disc area (Fabian et al 



2000). However, this may not be too surpris- 
ing since others (ISulentic et al.lll998al lbh have 



concluded that the evidence for and against 
an accretion disc origin for the Fe Ka emis- 
sion is equivocal. 

Chatteriee et al.1 f|201lh have reported that 



the fluctuations in the X-ray continuum and 
Fe Ka line are correlated. In their study these 
parameters were sampled approximately ev- 
ery 0.25 yr between 2004 and 2010. In order 
to examine the variations in these parameters 
more closely we have plotted them here in Fig 
7. To make the comparison easier the Fe Ka 
line values have been multiplied by a factor of 
0.6 to normalize the two curves. Both the (2-6 
Kev) continuum and the Fe Ka line intensity 
clearly show the same peaks. Furthermore, 
these are the same X-ray peaks found above to 
be correlated with the radio peaks and, with 
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a long delay, with the event start times. In 
each peak in Fig 7 the continuum intensity 
rises slightly before that of the Fe Ka line. 
A cross-correlation between these two lines is 



shown in Fig 8. As found bv lChatteriee et al 



( 20 111 ) this peak represents a significant cor- 
relation with a probability of P = 4 x 10 -4 . 
Fig 8 indicates that the Fe line lags the con- 
tinuum by ~ 50 days and it can be assumed 
that the variable portion of the line flux orig- 
inates in material located within ~ 50 lt-days 
of the (2-6 Kev) continuum radiation source. 

If the X-ray continuum radiation originates 
in the pc-scale jet, as found above, this evi- 
dence then indicates that the Fe Ka line radi- 
ation must originate nearby, and from the ex- 
cellent correlation found it seems likely that 
the Fe Ka line is somehow associated with 
the jet. It cannot originate inside the jet if 
the jet is moving relativistically, since the Fe 
line is not significantly shifted. It may origi- 
nate, however, from a bow wave produced by 
the ejected blob as it plows through an out- 
flowing wind. This might be expected to pro- 
duce some line broadening, but little line shift, 
as seen, since the speed of any outflow or wind 
will be quite low compared to that of the ma- 
terial in the jet. That the Fe Ka line radia- 
tion would be correlated with the appearance 
and disappearance of the continuum bursts as 
seen, would also be expected in this scenario. 

To the best of our knowledge there is no 
bursting radiation associated with the kpc- 
scale components that would have been signif- 
icant in the single-dish radio measurements. 
Therefore the bursting radiation must be as- 
sociated either with the accretion disc or the 
pc-scale jet. Because the X-ray telescope does 
not have the resolution required to resolve the 
accretion disc area from the pc-scale jet, this 
paper is the first to present observational evi- 
dence that can pinpoint the point of origin of 
the bursting X-ray component. 

Although the purpose of this study is to 



determine exactly where the bursting com- 
ponent of the X-ray radiation originates, we 
assume that the X-ray emission mechanism 
in the pc-scale jet is similar to that sug- 
gested for the external jets of powerful FR 
II sources which has been interpreted as in- 
verse Compton scattering off of cosmic mi- 
crowave background photons by the electrons 
in the relativistic jets. Recently, more and 
more evidence is being found showing that 
there is a connection between cm/mm ra- 
dio and high energy emission in quasars, es- 
pecially gamma-ray emission which has also 
been shown t o come from t he pc -scale jet. 
For example, ISchinzel et al.l (|201ll ) have re- 
cently shown that the gamma-ray emission in 
3C 345 is produced, not in a compact region 
near the central engine of the AGN, but in the 
Compton-loss dominated zone of the pc-scale 
jet. It should be noted that although these au- 
thors refer to 3C 345 as the archetypal quasar, 
this source has also been classified as an Ul- 



tra Lu minous Infrared Galaxy by iHou et al 
(|2009h . 



6. Conclusions 

Using a cross-correlation analysis we have 
found a strong correlation between the 14.5 
GHz radio and the X-ray light curves of 3C 
111, with the radio peaks trailing the X-ray 
by ~ 100 days. We have also demonstrated 
that this bursting component of the X-ray ra- 
diation originates, as does the bursting com- 
ponent of the radio radiation, in the ejected 
plasma blobs as they move outward in the pc- 
scale jet of 3C 111, peaking ~ 1 year after the 
start of each new event. The fact that a strong 
correlation exists between the start of each 
new radio event and its associated X-ray peak 
that occurs ~ 328 days later, indicates that 
the X-ray radiation must be associated with 
the radio blobs in the pc-scale jet that peak 
on this same time scale. Since the radio peak 
occurs ~ 100 days after the X-ray peak, it 
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will occur ~ 428 days (14.3 months) after the 
beginning of each new event. This is in good 
agreemen t with what was found for the radio 
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